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Abstract: A general procedure for construction of the formalism of quantum beam optics for any particle is reviewed. 
The quantum formalism of spin-^ particle beam optics is presented starting ab initio with the Dirac equation. As an 
example of application the case of normal magnetic quadrupole lens is discussed. In the classical limit the quantum 

q ■ formalism leads to the well-known Lie algebraic formalism of classical particle beam optics. 
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Whenever the possibility of a quantum formalism of particle beam optics is mentioned the immediate response, 
pH ' invariably, in the accelerator physics community is to ask what is the need to use quantum mechanics when classical 
O . mechanics has been so successful in the design and operation of numerous accelerators? Of course, this is a natural 
question and, though the system is quantum mechanical at the fundamental level, in most situations classical mechanics 
is quite adequate jjj since the de Broglie wavelength of the (high energy) beam particle is very small compared to the 
typical apertures of the cavities in accelerators as has been pointed out clearly by Chen. Q But, the recent attention to 
the sensitivity of tracking of particle trajectories to quantum granularities in the stochastic regions of phase space || 
and the limits placed by quantum mechanics on the achievable beam spot sizes in accelerators |I| clearly indicates the 
need for a formalism of quantum beam optics relevant to such issues. || Besides this, with ever increasing demand 
for higher energies and luminosity and lower emittance beams, and the need for polarized beams, the interest in the 
' studies on the various quantum aspects of beam physics is growing. || So, it is time that a quantum formalism of 
particle beam dynamics is developed in which all aspects (optical, spin, radiation, . . ., etc.) are considered in a unified 
framework. 



The grand success of the classical theories accounts for the very few quantum approaches to the charged-particle 
beam optics in the past. Notable among these are: 

• 1930 Glaser: Quantum theory of image formation in electron microscopy - Semiclassical theory based on the 



> 

in 
oo 
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^vq , nonrelativistic Schrodinger equation. J7[ 

• 1934 Rubinovicz; 1953 Durand; 1953 Phan-Van-Loc: Studies on electron diffraction based on the Dirac equa- 
tion. § 

• 1986 Ferwerda et ai: Justified the use of scalar (Klein-Gordon) equation for image formation in practical electron 
i ^ • microscopes operating even at relativistic energies. @ 

CZ2 ' 

• 1989-90 Jagannathan et ai: The first derivation of the focusing theory of electron lenses using the Dirac 
equation. B 1995: Quantum theory of aberrations to all orders using the Klein-Gordon theory and the Dirac 
Theory. || 1996: Spin dynamics of the Dirac particle beam. JTo| ] 



The formalism of quantum theory of charged-particle beam optics developed by Jagannathan et ai, based on the 
Klein-Gordan and Dirac equations, provides a recipe to work out the quantum maps for any particle optical system 



up to any desired order. JS| |11|. The classical limit (de Broglie wavelength — > 0) of this quantum formalism 



reproduces the well-known Lie algebraic approach of Dragt et al. [g_2| for handling the classical beam optics. Spin 
evolution, independent of orbital motion, can also be treated classically using the Lie algebraic approach. |l3| ] This 
brief note is to present the essential features of the quantum formalism of spin-i particle beam optics based on the 
Dirac equation. 



*To apear in the Proceedings of the 18th Advanced ICFA Beam Dynamics Workshop on Quantum Aspects of 
Beam Physics (QABP) (1 5-20 October 2000, Capri, ITALY) Editor: Pisin Chen (World Scientific, Singapore, 2002) 
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II. THE GENERAL FORMALISM OF QUANTUM BEAM OPTICS 



In many accelerator optical elements the electromagnetic fields are static or can be reasonably assumed to be static. 
In such devices one can further ignore the times of flights which may be negligible, or of no direct relevance, as the 
emphasis is more on the profiles of the trajectories. The idea is to analyze the evolution of the beam parameters of the 
various individual charged-particle beam optical elements (quadrupoles, bending magnets, • • •) along the optic axis of 
the system. Let us consider a charged-particle at the point (rj_, Si n ) where r± is the transverse coordinate and s refers 
to the coordinate along the optic axis. After passing through the system this particle arrives at the point (t"j_, s ou t)- 
Note that (r±, s) constitute a curvilinear coordinate system, adapted to the geometry of the system. Given the initial 
quantities at an Si n , the problem is to determine the final quantities at an s ou t, and to design an optical device in 
such a way that the relations between the initial and final quantities have the desired properties. Since we want to 
know the evolution of the beam parameters along the optic axis of the system the starting equation of the quantum 
formalism should be desirably of the form 

d 

ih—i(j(rx;s)=H^(rx;s), (1) 

linear in d/ds, irrespective of the basic time-dependent equation (Schrodinger, Klein-Gordon, Dirac, • • •) governing 
the system. So the step -I of building the quantum formalism is to cast the basic equation of quantum mechanics, 
relevant for the system under study, in the form ([!]) . Once this is done the step - II would be to obtain the relationship 
for the quantities at any point s to the quantities at the point Si n . This in the language of the quantum formalism 
would require to obtain the relationship for an observable {(O) (s)} at the transverse plane at s to the observable 
{(O) (sin)} at the transverse plane at Si n . This can be achieved by integrating ([[]). Formally, 

ip (r±;s) = U (s, s in ) ip (r±; s in ) , (2) 
which leads to the required transfer maps 

(O) ( Sin ) — > (O) (s) = (1> (s) \0\ ^ (s)) = (j> (,s in ) \&Ou\ V , ( Sin )^ . (3) 

Equation ([!]) is the basic equation of the quantum formalism of charged-particle beam optics and we call it as the 
beam optical equation, TL as the beam optical Hamiltonian and ip as the beam optical wavefunction. 

To summarize, we have a two-step algorithm to build a quantum formalism of charged-particle beam optics. On 
may question the applicability of the two-step algorithm: Does it always work? From experience we know that it 
works for the Schrodinger, Klein-Gordon and Dirac equations. The above description gives an oversimplified picture 
of the formalism than, it actually is. There are several crucial points to be noted to understand the success of the 
two-step algorithm. The first step in the algorithm to obtain the beam optical equation is much more than a mere 
mathematical transformation which eliminates V in preference to a variable V along the optic axis. There has to 
be a clever set of transformations ensuring that the resultant s-dependent equation has a very close physical and 
mathematical analogy with the original ^-dependent equation of the standard quantum mechanics. Without this 
guiding requirement it would not be possible to execute the second step of the algorithm which ensures that we can 
use all the rich machinery of the quantum mechanics to compute the transfer maps characterizing the optical system. 
This summarizes the recipe of obtaining the quantum prescriptions for the optical transfer maps. Rest is mostly a 
computational affair which is inbuilt in the powerful algebraic machinery of the algorithm. As in any computation, 
there are some reasonable assumptions and some possible approximations coming from physical considerations. It 
is important to note that in the case of the Schrodinger, Klein-Gordon and Dirac equations the beam optical forms 
obtained are exact. Approximations necessarily enter only in the step -II of the algorithm, i.e., while integrating the 
beam optical equation and computing the transfer maps for the quantum averages of the beam observables. As in 
the classical theory, the approximations arise due to the fact that only the first few terms are retained in the infinite 
series expansion of the beam optical Hamiltonian. The beam optical Hamiltonian is obtained as a power series in 
|7Tj_/po I where po is the design (or average) momentum of the beam particles moving predominantly along the optic 
axis of the system and 7Tj_ is the small transverse kinetic momentum. The leading order contribution gives rise to the 
paraxial or the ideal behavior and higher order contributions give rise to the nonlinear or aberrating behavior. Both 
the paraxial and the aberrating behaviors deviate from their classical nature by quantum contributions which are in 
powers of the de Broglie wavelength of the beam particle (Ao = 2ttTi/pq). The classical formalism is obtained from 
the quantum formalism by taking the limit Aq — ► 0. 
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III. FORMALISM OF THE DIRAC PARTICLE BEAM OPTICS 



Now we shall see how the above algorithm works for the Dirac particle. Let us consider a monoenergetic beam of 
Dirac particles particles of mass m, charge q and anomalous magnetic moment transported through a magnetic 
optical element with a straight optic axis characterized by the static potentials (4>(r), A(r)). The beam propagation 
is governed by the stationary Dirac equation 

&D\rl> D )=E\il> D ) , (4) 

where \iPd) is the time-independent 4-component Dirac spinor, E is the total energy of the beam particle and the 
Hamiltonian H^, including the Pauli term, is given by 

= (3mc 2 + ca ■ (-ihV - qA) - /i a/ 3£ B , (5) 

where the symbols have their usual meanings. [Q To cast ([|) in the required beam optical form ([!]) we multiply Hd 
(on the left) by a z /c and rearrange the terms to get 

d 

H D = -poPx^z - qA z I + a z otj_ ■ ttj_ + (/x / 'c)/3a 2 £ B , (6) 



where x is a diagonal matrix with elements — I/O — V£) an( l £ = Vv^ + mc 2 )/(E — mc 2 ). Equation (||) is still 
not in a completely desirable form. So we resort to a further transformation: 

— N0 = MhM, M=^=(I + X a z ). (7) 

Then we obtain 

ih 4- W) = H' W) , H' = MHdM- 1 = - Po /3 + S + 6, (8) 
oz 

where the nonvanishing matrix elements of the even term £ and the odd term O are given by 

£ n = -qA z t - ( Pa /2c) {(£ + r 1 ) <Tx-B x + (£- C 1 ) , 

= - (m /2c) { (e + r 1 ) <r±-B±-(€- r 1 ) , 

012 = f [o- J_ • 7T± - (A*o/2c) {l (f - (-B a O-„ - Byd,) 

-(e + r 1 )^!}] , 

021 = -r 1 [tr ± ■ tt ± + ( Ma /2c) {« (£ - r 1 ) {B x o y - B y a x ) 

+ (t + C 1 )B z t}] . (9) 

The effect of the transformation (^) is to make the lower components of a Dirac spinor corresponding to a quasi- 
paraxial beam moving in the positive z-direction negligible compared to the upper components and thus effectively 
making the 4-component spinor as a 2-component spinor. Now one may observe the close analogy: 

Standard Dirac equation Beam optical form 

mc 2 f3 + £- D + 6 D -poP + S + 6 

Positive energy Forward propagation 

Nonrelativistic, c |7r| <C mc 2 Paraxial beam, |7Tj_| <C po 

mc 2 : Note ih^- w mc 2 Tp -p : Note i?i§§ m -p ip 

Nonrelativistic motion Paraxial behavior 

+ Relativistic corrections + Aberration corrections 

This completes the step -I of the algorithm. To execute the step -II we proceed as follows. 

The above analogy suggests that, as the most systematic way to understand the Dirac Hamiltonian as a nonrel- 
ativistic part plus relativistic correction terms is to use the Foldy-Wouthuysen (FW) transformation technique, |l4| ] 
we should adopt a similar FW-like approach in the beam optical case to understand the beam optical Hamiltonian 
as a paraxial part plus nonparaxial correction terms. This leads to a procedure to obtain the paraxial behavior 
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accompanied by a systematic method to compute the aberrations to all orders in powers of the expansion parameter 
1/po- To leading order, the first FW-like transformation is 



(i) 



Then 



az 



1> 



1 



exp(-/3a/2p W 



6<n 



(10) 



Po (3 + £ (1) +6 {1) , 
d 



ih—o 

oz 



(11) 



£(i) = g_ J_po^ + .. OW = -— ^ 0,£ 
2p 2p Li- 

lt is to be noted that the transformation (jl^) keeps the upper components of the beam optical wavefunction large 
compared to its lower components. One can proceed with further FW-like transformations and stop at any desired 
stage. Let us denote the 2-component spinor comprising the upper components of the final 4-component spinor 
obtained in the above process as \ip). 

Up to now, all the observables, the field components, time etc., have been defined in the laboratory frame. The 
covariant description of the spin of the Dirac particle has the simplest operator representation in the rest frame of 
the particle. Thus, in accelerator physics the spin is defined in the rest frame of the particle. So we make a further 
transformation which takes us from the beam optical form to the accelerator optical form 



4> {A) \ = exp 



1> 



(12) 



Thus, up to the paraxial approximation the accelerator optical Hamiltonian Jl^,|ll| is 

d 



ih- 



dz 



H(A) 



-po - qA z + —ttI + -^Q s ■ S , 



2 Po 



Po 



(13) 



where O s = — ^ {qB + e (B\\ + jB ± )}, 7 = E/mc 2 and e = 2m/i a /?i. 



IV. AN EXAMPLE OF APPLICATION: MAGNETIC QUADRUPOLE LENS 



Let an ideal normal magnetic quadrupole of length £, characterized by the field B = (—Gy, —Gx, 0), be situated 
between the transverse planes at z = z- ln and z = z out = z- m + £. The associated vector potential can be taken to be 
A= (0,0,iG(x 2 -2/ : 
is 



! )) with G as constant inside the lens and zero outside. The accelerator optical Hamiltonian |l(J 



'P0 



H(z) 



1 Pi 



•2po 



H L {z) = -pa + -±-pl - \qG (x 



for 



2p ^-L 2 



for z < z; n and z > z ont , 

V 2 ) +*¥ L (V<T* + XVy) 

with rj = (q + je)G£h/2pl . 



(14) 



The subscripts F and L indicate the field-free and the lens regions respectively. 

Best way to compute the ^-evolution operator U is via the interaction picture, used in the Lie algebraic formu- 
lation |L2| of classical beam optics. Using the transfer operator thus derived |ll| we get the transfer maps for the 
averages of the transverse phase-space components: with the subscripts in and out standing for (z in ) and (z out ), 
respectively, 
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Tn = M- > MoM < 



Mr 



( cosh(Vl«) 
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-A=sinh(v^) 
cosh(y/K£) 
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(15) 
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1 J 






cos(V^) 
Ksm{yfKt) 






^sin(^) 
cos(VK£) / 



Thus we have got a fully quantum mechanical derivation of the combined effect of the focusing action of the quadrupole 
lens (note the traditional transfer matrices) and the Stern-Gerlach force. It may be noted that the quantum formalism 
of spinor beam optics supports, in principle, the idea of a Stern-Gerlach spin-splitter device to produce polarized 
beams. ]l5[ ] The transfer map across the quadrupole lens for the spin components computed using the above accelerator 
optical Hamiltonian describes the well known Thomas-Bargmann-Michel-Telegdi spin evolution, pffl 



V. CONCLUDING REMARKS 



In fine, we have seen how one can obtain the formalism of quantum beam optics for any particle, starting ab initio 
from the relevant basic quantum equation, at the single-particle level. A two-step algorithm for this purpose has been 
suggested. Using the general principle, the construction of a spinor theory of accelerator optics, starting from the 
Dirac equation and taking into account the anomalous magnetic moment, has been demonstrated. As an example of 
application of the resulting formalism the normal magnetic quadrupole lens has been discussed. In the classical limit 
the quantum formalism leads to the Lie algebraic formalism of charged-particle beam optics. 

To get a formalism taking into account the multiparticle effects, particularly for the intense beams, it should be 
worthwhile to be guided by the quantum- like approaches to the particle beam transport: Thermal Wave Model |H| 
and Stochastic Collective Dynamical Model p7| . Recently the quantum- like approach has been applied to construct a 
Diffraction Model for the beam halo. [fl8f This model provides numerical estimates for the beam losses. In this context, 
another useful approach could be to use the Wigner phase-space distribution functions. Hcinemann and Barber |l9| ] 
have initiated the derivation of such a formalism for the Dirac particle beam physics starting from the original work 
of Derbenev and Kondratenko |2(| who used the FW technique to get their Hamiltonian for radiation calculations. 

The present study is confined to systems with straight optic axis. An extension to the curved optic axis systems 
should be done. This would involve the subtlities of quantization in curvilinear coordinates. Then there are the well 
known questions related to the position operator in the relativistic quantum theory. Also, there are doubts about the 
exact form of the Stern-Gerlach force for a relativistic particle. pl| To address such questions from the point of view 
of experiments using particle beams the right platform would be the formalism of quantum beam optics. 
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